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We study the perturbation expansion of the free energy of A/" = 4 supersymmetric 
SU{N) Yang-Mihs at finite temperature in powers of 't Hooft's couphng g'^N in the 
large N hmit. Infrared divergences are controUed by constructing a hierarchy of two 3 
dimensional effective field theories. This procedure is applied to the calculation of the 
free energy to order (g'^N)^''^, but it can be extended to higher order corrections. 
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In recent years some light has been shed on understanding the relation between string 
theory and gauge theories |]l|. A particularly interesting system is A^ coincident, parallel 
D3-branes. This system realizes N" = 4 supersymmetric U{N) Yang-Mills in its world vol- 
ume (4 dimensions). In the large N limit, the 3-brane system becomes a black brane whose 
Bekenstein-Hawking entropy can be obtained by considering Ramond-Ramond charged 3- 
brane classical solutions . It is therefore interesting to compare the 3-brane thermodynam- 
ics with that of A/" = 4 supersymmetric U{N) Yang-Mills. The free energy density of A/" = 4 
supersymmetric U{N) Yang-Mills in the ideal gas approximation is ^ideai = —ii'^N'^T^/Q. On 
the other hand, the free energy of a black 3-brane was found to be J^br = (3/4) x J-'idea.i- 
Maldacena's conjecture helps to understand the relative factor 3/4 between J-'ideai and 
J^BR- The conjecture relates type IIB superstrings on AdSs x S^ to A/" = 4 supersymmetric 
SU{N) Yang-Mills in 4 dimensions and allows to compute the next-to-leading order correc- 
tion to J^BR in the strong coupling limit @], 
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where A = g'^N is the 't Hooft coupling. /(A) = J^/J-'ideai should be interpreted as a function 
whose strong coupling limit is f{oo) = 3/4, while the weak limit is /(O) = 1. It has been 
suggested that /(A) is a monotonic function that interpolates between the strong and the 
weak coupling limits |]I], Q. 

We will explore the weak coupling limit of A/" = 4 supersymmetric SU{N) Yang-Mills 
(SYM for short) in 4 dimensions at finite temperature (3 space dimensions). In the present 
note, we will carry out the calculation of /(A) up to order A^'^ by using the effective field 
theory approach at finite temperature |^ |]. Related work can be found in and @, 0, where 
the terms of order A and A^/^ have also been computed, respectively, but using approaches 
different from ours. 

Feynman rules at finite T are the same as those at T = except that loops involve infinite 
sums over Fourier modes |Tn]. The thermodynamic properties of A/" = 4 supersymmetric 



SU{N) Yang-Mills in 4 dimensions, in thermal equilibrium are described by the free energy 
density JF = —T\ogZ/V, where Z is the partition function and V is the 3d volume of the 
system, which is taken to infinity in the thermodynamic limit. The partition function can 
be written as a path integral over the fields of the Lagrangian C that describes SYM. Fields 
are in the adjoint representation of the gauge group, the gauge coupling is g. An expansion 
for — T log Z/V in powers of the coupling constant g can be obtained by summing up the 
vacuum Feynman diagrams that contribute at a given order in g. 

The infrared (IR) behavior of field theories at finite T however prevents us from computing 
the free energy density in the form just described. There are two sources of IR divergences. 
First, some fields that are massless at T = 0, become massive at T 7^ with a mass of 
order gT. Diagrams with self-energy insertions have IR divergences that do not cancel order 
by order in the coupling constant. A second source of IR divergences, which is specific to 
non-abelian gauge theories, has to do with the gauge field self-interaction terms. 

The effective field theory approach allows to remove systematically the IR divergences 
associated to thermal mass insertions. Also, IR divergences related with the gauge field 



self-coupling terms are isolated and interpreted as a source of non-perturbative effects [11 . 



In the effective field theory approach for gauge theories, a hierarchy of two effective field 
theories in 3 dimensions which only contain static modes are constructed. Since, in addition, 
these effective field theories do not contain fermionic fields, they are very convenient for 
non-perturbative studies on the lattice [§].0 

If non-static modes are integrated out in the partition function Z^ we can write 

Z = J VAo{^)VA{^)VM^) e-/'^'^^°« , (2) 

where £efr is an effective Lagrangian compatible with the internal symmetries of C. C^s 
describes a 3 dimensional effective field theory |T^. £eff is made out of an electrostatic field 



Aq[x), a magnetostatic gauge field v4j(x), and 6 scalars 0/(x), where / = 1, ... ,6. These 
fields can be identified, up to field redefinitions, with the static modes of their 4 dimensional 
counterparts. Note that fermion fields have been integrated out completely because they do 
not have a static mode. We can write £eff = /e + '^esym, where /e is a constant that cannot 
be ignored when computing the free energy. C-esym contains the remaining terms and has 
the form 

/^ESYM = \G1^Gt^ + ^{D.AonD.AoT + ^ (A0/)'^(A0/)'^ 

+ ^mlA^AS + ^mlriri + SCesym, (3) 

where G^j = diA'^—djA°^+gEf°''''^A'lA'^ is the magnetostatic field strength with gauge coupling 
constant qe- The term SCesym in (§|) represents an infinite number of terms constructed 
out of Aq{-x), Aj(x), and 0/(x) that contribute only at order higher than A^/^. SCesym 
contains non-renormalizable interaction terms. We call the effective theory constructed by 
integrating out non-static modes Electrostatic SYM (ESYM). The free energy can be written 
as J^ = T{fE -log Zesym/V), where 

^ESYM = J VAo{^)VA{^)VM^) e-Z'^'^^-s™ . (4) 

We can go a step further by integrating out Ao(x), and </)/(x) in @. We obtain 

2esym= fvAi^)e'I'^'^^, (5) 



where Ces is an effective Lagrangian compatible with the internal symmetries of the theory 
which is made out of a magnetostatic gauge field Aj(x). We can write C^s = Im + >Cmsymj 
where Jm is a constant. £msym is of the form 

-CmSYM = T ^ij^ij + "^-^msym, (6) 



■^Several theories at finite temperature have been studied perturbatively and non-perturbatively using the 
effective field theory approach. The ^"'-theory is studied in g, QCD in jl^, |ll|, QED in [l3| , and scalar 
QED in [|lj|. The electroweak phase transition has been studied in the Standard Model |^, |l5|] and in some 
of the Minimal Supersymmetric extensions of the Standard Model p6| . 



where Hij = d^A"^ — djA°- + gMf"''''^A\A^- is the magnetostatic field strength with couphng 
constant Qm- The term 5£msym in (U) represents an infinite number of terms constructed 
out of Aj(x). We call this theory Magnetostatic SYM (MSYM) because it is made out of 
magnetostatic fields only. The free energy can therefore be written as 

^ = T(/. + /M-i^^^^), (7) 

where 

^MSYM = JVA,{^) e-/<^^-^MSYM . (8) 

We may think of using (J^) to compute the free energy density of SYM. We need to 
determine Z^;, the effective parameters of £esym, /a/, and the effective parameters of £msym- 
The term Je and the effective parameters of £esym may be determined by computing the 
same physical quantities in both the full theory (SYM) and in the effective field theory 
described by £eff (ESYM) and comparing the result. At first sight this procedure may seem 
useless if we consider that to compute physical quantities in the full theory we have to remove 
the IR divergences. It is like having to face the very same problem we wanted to solve at the 
beginning. However, to compare results, it is actually enough to compute them in a region 
where we know that both theories describe the same physics. 

As mentioned above, self-energy insertions of order gT in the full theory give rise to IR 
divergences. If we introduce an IR momentum cutoff A^;, we will obtain physical quantities 
as expansions in powers of gT j Ke- The expansion diverges when A^; goes to zero. However, 
if we demand A^; ^ gT ^ the expansion makes sense. Of course, the result depends on A^; and 
one cannot get rid of it unless an infinite number of diagrams are summed up. In the effective 
theory, the thermal masses {vhe and rus) are also of order gT (we will explicitly see it below). 
Therefore, an IR momentum cutoff A^; ^ gT is equivalent to treating the mass terms in 
£eflF as interaction terms. By computing physical quantities in this way and comparing the 
results, we are able to determine the parameters of the effective theory (ESYM) as functions 
of T, (yf, and A^; 3> gT . 

Once we have determined the parameters of ESYM by matching, we can use ESYM as an 
ordinary field theory. In general, when computing diagrams using ESYM, we will encounter 
ultraviolet (UV) divergences. These UV divergences can be regulated by introducing a UV 
momentum cutoff A. Moreover, ESYM is not fully free of IR divergences. In order to 
regulate the IR divergences originated by the gauge boson self-interaction terms in ESYM, 
we have to introduce an IR momentum cutoff Km- Analogously to what happened to the 
full theory with IR cutoff A^;, we can make sense of the perturbative expansion of ESYM 



only if Am > g^T |Tg, ^. Note that the UV cutoff of ESYM has to satisfy A > Am- Since 
Ag 3> gT and Am ^ g'^T, we can think of A^; as A^; ~ T and Am as Am ~ gT. Therefore, 
we can take A^; as the ultraviolet cutoff of ESYM: A = A^;. The A^ dependence of the 
parameters is canceled by the A^ dependence of the loop integrals in ESYM. 

We see from (^ that MSYM is a pure gauge theory in 3 dimensions. It is a confining 
theory that cannot be studied perturbatively. Its Feynman diagrams only contain gauge field 
self-interaction terms and one should expect the same IR problems that we have already 
pointed out for the full theory and ESYM. However, if we introduce an IR momentum cutoff 



^&i ^ g'^T, we can still use a diagrammatic expansion to determine the effective parameters 
of MSYM by matching physical quantities with ESYM. Of course, although we may be 
able to determine the effective parameters of MSYM, — log^MSYM/^ in (0) has still to be 
calculated non-perturbatively. 

We write "X ^ ■ ■ ■" to denote that the physical quantity X has been computed using 
the IR cutoffs A^; or Am- An expression of this type for X is suitable for matching but it 
should not be confused with a correct perturbative expansion (which should have the IR 
divergences removed). The perturbative expansions performed with IR cutoffs are called 
strict perturbation expansions. 

So far, we have used momentum cutoffs to learn how to organize the perturbative ex- 
pansion of A/" = 4 supersymmetric Yang-Mills at finite T. However, it turns out to be more 
convenient to use a regulator that automatically gets rid of power divergences. Such a reg- 
ulator is dimensional regularization. We therefore consider the full theory in (3 — 2e) -|- 1 
dimensions. Since, in the present note, we are interested in computing the free energy den- 
sity of A/" = 4 supersymmetric SU{N) Yang-Mills to order A^/^, we limit ourselves to the 
effective parameters and corrections that contribute to such an order. 

The strict perturbation expansions in the full theory can be greatly simplified if we 
consider that A/" = 4 supersymmetric Yang-Mills in 4 dimensions (SYM4) can be obtained by 
dimensional reduction of A/" = 1 supersymmetric Yang-Mills in 10 dimensions |19[ . Physical 
quantities in SYM4 can be obtained perturbatively by using SYMio Feynman rules while 
loop integrals are performed in 4 dimensions (SYMio^4) @]. In this formulation, the lOd 
Afj,{fi = 0,1, 2, 3) can be identified, up to scale redefinitions, with the 4d gauge field A^. The 
remaining A/+3 (/ = 1, . . . , 6) are identified with the 6 scalars </>/ of A/" = 4 supersymmetric 
Yang- Mills in 4d. While leading order contributions to the strict perturbation expansions in 
the full theory can be easily computed using either SYM4 or SYMio^4, higher than leading 
order contributions are much less arduous in SYMio^4. 

The gauge coupling constant of ESYM, qe, can be read off from the Lagrangian of the 
full theory. By substituting Aq{x, r) -^ \/TAq(x) in the Lagrangian of SYM4 and comparing 
/q drCsYM with /^esym, we find out, to leading order in g'^, 

gl = g'T . (9) 

In QCD, the gauge invariant electric screening mass is given by the pole location of 
the /i = z/ = component of the complete gluon propagator [^. We will use a sim- 
ilar definition in SYM. Using SYMio^4, the inverse gauge field propagator has the form 
(^"^ [Sfiuk'^ + n^i,(0, k)]. In order to compute the static masses in £esym, it is convenient to de- 
fine nei(fc2) and nsc(A;2) such that noo(0, k) = U,i{k^) and n,+3,j+3(0, k) = Sijl^scik^) (/, J = 
1, . . . , 6), respectively. Calculations are performed in Feynman gauge. 

The electrostatic mass tue is determined by computing the electric screening mass in 
the full theory and matching it with the result obtained in ESYM. The electric screening in 
SYM is the solution to the equation 

k'^ + U^i{k^) = atk'^ = -ml. (10) 
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Figure 1: (a) One- loop Feynman diagrams for the gauge field self-energy. (6) Two-loop 
Feynman diagrams for the free energy. Wavy lines, solid lines, and dotted lines represent 
the propagators of gauge fields, fermions, and ghosts, respectively using SYMio^4. 

Note that (|T^) gives rise to a double expansion in powers of g (or equivalently, in number of 
loops) and in powers of m^i ~ g^T"^- In the calculation of the free energy to order A^^^, it is 
enough to compute Ilei at one-loop order. Then, the expression for the electric screening mass 
to leading order in g"^ is m1^ Ri Ilgj (0). Diagrammatically, Ilgj (0) is given by the graphs 
shown in Fig. |l|(a). These diagrams give rise to sum-integrals that can easily be evaluated 



by using standard methods [T^. We find that, in the full theory, the strict perturbation 
expansion for rriei is 

ml^2CAg^T\ (11) 

In ESYM, the electric screening mass rriei gives the location of the pole in the propagator 
for the field AqIx). Denoting the self-energy function by IlE{k'^)S°'^, "n^ei is the solution to 

A;^ + m| + UEik"^) = at A;^ = -m^. (12) 

Similarly to (pi]|), (|T^) also gives rise to a double expansion in number of loops and in 
powers of mli ~ g^T"^- The strict perturbation expansion of rriE is obtained as an expansion 
involving Feynman diagrams evaluated at zero external momentum. Since ttie, rns, and 
gE are treated as perturbation parameters, there is no energy scale in the dimensionally 
regularized integrals and they all vanish. The solution to the equation ([12D for the screening 
mass is therefore trivial: m^i ~ m"^. Comparing this result with ([TT|), we find that, to 
leading order in g"^ 

ml = 2CAg'T\ (13) 

The calculation of the scalar static mass ms is similar to the calculation of tue- We 
define a scalar screening mass rUsc which satisfies fc^ + Uscik"^) = at A;^ = —m'^c i^ the full 
theory and an equation analogous to ([T^ ) in the effective theory. After comparing the result 
in SYM with the result in ESYM, we find that, to leading order in g"^ 

ml = CAg'T' . (14) 

The effective parameter Je is evaluated by matching the calculation of the free energy 
density in the full theory and in the effective theory. The leading order contribution to the 



free energy density in the full theory is given by a familiar result from blackbody radiation 
j^W ~ —TT^T^^gg ^^^ _|_ y^^/g^^ where 5b {5p) is the number of bosonic (fermionic) degrees 
of freedom of the theory. In the case of A/" = 4 SYM4, 5b = 5f = Sc?^, where (Ia is the 
dimension of the adjoint representation of SU{N): cIa = N"^ — 1- Alternatively, J-"*^^-* may 
be obtained by computing the one-loop vacuum diagrams of SYM4 using SYMio^4. The 
next-to-leading order contribution to the free energy density is given by the 2-loop Feynman 
diagrams shown in Fig. |I](6) using SYM10-+4. The evaluation of these diagrams involves 



elementary sum-integrals that can easily be evaluated by standard methods |T^. Adding 
the 2-loop order contribution to JF^^), we obtain that the strict perturbation expansion of 
the free energy density in the full theory is 



t^'^cLa^a /. 3 



r^-^T-[i--c.,'y (15) 

In ESYM the effective parameters irtE-, m,s, and qe sue treated as perturbation param- 
eters. Then, as there is no scale in the dimensionally regularized integrals, they all vanish. 
The strict perturbation expansion of the free energy density is therefore trivial: J-' ~ TJe- 
Comparing this result with (|15D , Je to order g"^ is 

As a check we have also computed the effective parameters of ESYM ttle, ttls, and J'e using 
SYM4 instead of SYMio^4 and found the same results given in Eqs. (|T3|), (^^, and (|I^) 



To determine /m, we have to compute the strict perturbation expansion of — log Z-esym/V 
and compare it with that of /m — log-^MSYM/^- Computing the one-loop vacuum diagrams 
of ESYM, we find — log^ESYM/^ ~ — c^a/IIStt) (m|; -|- 6 m|). The effective parameters 
of £msym are treated as perturbation parameters. Feynman diagrams give rise to integrals 
which have no scale dependence. Therefore, using dimensional regularization all the con- 
tributions to — log^MSYivi/V^ vanish. We have /m — ^ogZusYu/V ~ /m- Comparing this 
result with —logZESYu/V, we obtain /m'- 

/M = -^(m| + 6m|). (17) 

The gauge coupling constant qm can be read off from the Lagrangian of ESYM £esym- To 
leading order in g% and g"^, we obtain gj^ = 9% = 9^T. 

The free energy of A/" = 4 supersymmetric SU{N) Yang-Mills is given by (^. The term 
— log^MSYM/^ is non-perturbative. It is easy to find out at which order non-perturbative 
corrections become relevant. The only quantity with dimensions in £msym is gM (neglecting 
(^-^MSYM whose contribution is subleading). Since — logZusYu/V has dimension 3 in units of 
energy, its lowest order contribution to the free energy is therefore of order gf^ = (g'^T)^. We 
conclude that the non-perturbative corrections generated by log^MSYivi/^ only contribute 
to the free energy of SYM at order A^ or higher. 

Using (1^), the free energy of A/" = 4 supersymmetric SU{N) Yang-Mills in 4 dimensions 
to order A^/^ is given by J^ = T [/e - d^/(127r) (m| + 6m|)]. Substituting (|T3|), (0), 

6 



and (|T6|), we finally find that, in the large A^ limit, 

/(A) = /- = 1-^A + ^A3/^ (18) 

This result is in agreement with previously reported calculations. In Ref. 0, the free 
energy was evaluated to order A. In p[ and 0, the term of order A^/^ has been computed 
by regulating the IR divergences of SYM4 by adding and subtracting out thermal masses to 
the Lagrangian. Vazquez- Mozo used SYMio^4 to compute the first two terms in (p^. 

Since the term of order A^^^ in ( |T8| ) is positive, it seems that our result does not favor a 
monotonically decreasing function interpolating between the weak and strong coupling limits. 
However, such a behavior may be an effect of the zero convergence radius of the weak-coupling 
expansion (|TH]). It has been suggested [^ that there could be a phase transition in A at 
large A^. Even though perturbative calculations cannot settle the debate on this matter, it 
may be possible to find an indication of a phase transition by using Pade approximants. An 
analysis of (|1^) using Pade approximants has been carried out in 0. There, an indication 
has been found of a smooth interpolation between the weak and strong coupling regimes. 

The effective field theory approach provides a systematic way to compute perturbative 
contributions to the free energy of SYM while taking advantage of dimensional reduction 
(SYMio^4) to determine the effective parameters of ESYM. This is a major simplification for 
computing /(A) to order less than A'^ ||2^. Starting at order A'^, there are non-perturbative 
contributions as well as perturbative contributions. Non-perturbative effects are described 
by MSYM, a pure gauge theory in 3 dimensions. As we have already pointed out, the 
effective field theory approach gives rise to effective field theories in 3 dimensions, without 
fermions, that are suitable for non-perturbative studies on the lattice. Lattice simulations 
of the effective field theories could therefore be used to explore the non-perturbative sector 
of A/" = 4 supersymmetric SU{N) Yang-Mills. 
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